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ABSTRACT 

The pulsar PSR B1259— 63 is in a highly eccentric 3.4-yr orbit with the Be star SS 
2883. Timing observations of this pulsar, made over a 7-yr period using the Parkes 
64-m radio telescope, cover two periastron passages, in 1990 August and 1994 January. 
The timing data cannot be fitted by the normal pulsar and Keplerian binary parame- 
ters. A timing solution including a (non-precessing) Keplerian orbit and timing noise 
(represented as a polynomial of fifth order in time) provide a satisfactory fit to the 
data. However, because the Be star probably has a significant quadrupole moment, we 
prefer to interpret the data by a combination of timing noise, dominated by a cubic 
phase term, and u> and x terms. We show that the ui and x are likely to be a result 
of a precessing orbit caused by the quadrupole moment of the tilted companion star. 
We further rule out a number of possible physical effects which could contribute to 
the timing data of PSR B1259— 63 on a measurable level. 
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1 INTRODUCTION 

The pulsar PSR B1259— 63 is a member of a unique binary 
system. Discovered using the Par kes telescope in a sur vey 
of the Galactic plane at 1.5 GHz (Johnston et al. 1992a), it 
was shown by Johnston et al. (1992b) to be in a highly ec- 
centric 3.4-yr orbit with a lOth-magnitude Be star, SS 2883. 
The pulsar period, P, is relatively short, 47.8 ms, and the 
measured period derivative, P, gives a pulsar characteristic 
age, t c = P/(2P), of 3.3 x 10 5 yr and a surface magnetic 
field of 3.3 x 10 11 G. This therefore is a young system, which 
may evolve through an accretion phase to form a single or 
binary millisecond pulsar. Rapidly spinning neutron stars 
can only accrete matter if the co-rotation velocity at the 
Alfven radius is less than the Keplerian velocit y at the same 
radius (Bhattacharya & van den Heuvel 1991). Equality of 



these velocities defines the 'spin-up line'. At present, PSR 
B1259— 63 lies well to the left of the spin-up line, so that 
accretion onto the neutron star is not possible until either 
the pulsar slows down or the pulsar magnetic field decays. 
Timing observations of PSR B1259— 63, made over a 



3.4-yr interval and covering the 1990 August periastron, 
were reported by Johnston et al. (L994). A phase-connected 
fit to these data gave parameters for the pulsar and its or- 
bit, and showed that the next periastron would occur on 
1994 January 9. This paper also reported optical observa- 
tions which indicate that the companion star is of spectral 
type B2e, with a mass of m* ~ 10 Mq and radius J?* ~ 6 
Rq. The break-up velocity at the equator, iWx, for B2e 
stars is not very well known; it is estimated to be ~ 380 km 
s" 1 by Slettebak et al. (1980) and ~ 480 km s" 1 by Schmidt- 
Kaler (1982). Recent work by Porter (1996) suggests that 
most Be stars rotate at ~ 70 per cent of the break-up veloc- 
ity. From the mass function, a companion mass of 10 Mq 
and a pulsar mass, m p , of 1.4 Mq imply an orbital incli- 
nation i = 36°. The orbital eccentricity is very high, 0.87, 
and the pulsar approaches within 24 R* of the companion 
star at periastron, passing through the circumstellar disk. 
Extensive observations of the pulsar were made at several 
radio frequencies before and after the 1994 January perias- 
tron, in order to probe the circumstellar environment of SS 
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2883 (Johnston et al. 1996; Melatos, Johnston & Melrose 
1995). 

Observations made between 1990 January and 1994 Oc- 
tober were well explained by step changes in the pulsa r pe- 
riod at the two periastrons (Manchester et al. 1995), at- 
tributed to a propeller-torque spin-down caused by the in- 
teraction of the pulsar with the circumstellar matter at the 
Alfven radius (Illarionov & Sunyaev 1975, King & Cominsky 
1994, Ghosh 1995). 

In this paper we report on additional timing observa- 
tions made using the Parkes radio telescope over the past 
two years which, together with the earlier data, give a total 
timing data span of seven years. We find that the timing 
solution of Manchester et al. (1995) does not fit the recent 
data and discuss alternative models and their interpretation. 



2 TIMING OBSERVATIONS AND DATA 
ANALYSIS 

A total of ~300 pulse times of arrival (TOAs) were mea- 
sured at the Parkes radio telescope between 1990 January 
and 1996 December. Most of the observations were at fre- 
quencies around 1.5 GHz, but observations at 0.43, 0.66, 4.8 
and 8.4 GHz were also made. At all frequencies, dual-channel 
cryogenically cooled systems receiving orthogonal linear po- 
larizations were used. After conversion to an intermediate 
frequency, signals for each polarization were split into sub- 
bands using filterbanks with channel widths of 0.125 or 0.25 
MHz for frequencies below 1 GHz, and 5 MHz for higher fre- 
quencies. Signals from corresponding filter polarization pairs 
were detected, summed, high-pass filtered and one-bit dig- 
itized, usually with a sampling interval of 0.6 ms. Further 
detail s of the observing systems are given in Johnston et al. 
( |1996[ ). 

The data were folded at the topocentric period to form 
mean pulse profiles. Integration times were typically 10 min 
at 1.5 GHz, and 20 to 30 min at other frequencies. Data 
taken prior to 1994 were convolved with a 'standard profile' 
in the time domain to produce TOAs whereas data taken 
after 1994 were processed using a different method which 
involves matching the Fourier components of the observed 
and standard profiles. Uncertainties in the TOAs are some- 
what smaller using the second method. However in order not 
to bias the fitting procedure towards the second half of the 
data (i.e. the second orbit), we use unweighted fits through- 
out this paper; that is, all TOAs were assigned the same 
weight regardless of the estimated TOA uncertainty. Also, 
we use only TOAs obtained at frequencies around 1.5 GHz 
and those TOAs at 4.8 GHz, which show the best signal-to- 
noise ratio, in order to obtain a good value for the dispersion 
measure. Pulsar and binary parameters wer e obtained using 



the least- squares fitting program TEMPO (Taylor & Weis- 



berg 1989) with the Jet Propulsion Laboratory solar-system 



ephemeris DE200 (Standish 1982). TEMPO was extended 



by the addition of a new timing model for binary pulsars 
that orbit a companion star with a significant quadrupole 
moment (Wex 1998). 

As discussed by Johnston et al. fll996|) , significant dis- 
persion and scattering changes were observed around perias- 
tron. Because of this, data from 1990 July and 1993 Decem- 
ber were omitted from the analysis. The pulsar was eclipsed 



o ° 

3 
T3 



! 1 1 ' 

• 


1 1 | 1 1 1 1 | 1 1 


* i i i * * i i * * i i . 

! ' 


i 

( 

7,1. 





* '. 

■ i i ■ 



-2000 



-1500 



-1000 



-500 



500 



MJD-50000.0 



Figure 1. Post-fit residuals from an unweighted fit for pulsar 
position, period, period derivative, dispersion measure and Kep- 
lerian orbital parameters. The rms residual is 2100 fis. The ver- 
tical dotted lines in this and subsequent figures indicate the two 
epochs of periastron passage. 



during 1994 January. Thus, the two largest gaps in the data 
in terms of true anomaly, tp, are one of 129 days around the 
first periastron (tp = —145° . . . 100°) and 76 days around the 
second periastron (tp = —128° . . . 105°). This will emerge to 
be a problem when searching for the correct timing model 
for PSR B1259-63. 

The data were first fitted for pulsar position, period, pe- 
riod derivative, dispersion measure, and the five Keplerian 
orbital parameters. The best results give an rms residual of 
2100 /is and the post-fit residuals shown in Fig. 1. System- 
atic variations in the residuals are observed at all orbital 
phases, showing that this set of parameters does not model 
the timing behaviour of the system satisfactorily. 

Three types of model are considered in the following 
subsections with a view to minimising these residuals and 
understand their origin. 

2.1 Period steps at periastron 

As mentioned in the Introduction, observations made be- 
tween 1990 January and 1994 October were well fitted by 
step increases in the pulsar period at the two periastrons. 
Fig. 2 shows the pre-fit residuals ob tained using the timing 
solution given in Manchester et al. (1995) on the extended 



data set. Clearly, the timing observations made in 1995 and 
1996 are not modelled by this solution. 

There is other evidence which suggests that propeller 
spin-down does not occur at periastron. Models of the sys- 
tem by Tavani & Arons (1997) based on extensive X-ray and 
7-ray observations during the 1994 periastron (Kaspi et al. 
1995, Hirayama et al. 1996, Grove et al. 1995) conclude that 
accretion or propeller-powered processes are ruled out. 

If we allow fitting for a period step at the two pe- 
riastrons, we find a rather good fit, with rms residual of 
340/is. However, the two period steps have opposite sign, 
Av = -1.7xl0~ 8 at the first periastron and Av = 2.6xl0~ 8 
at the second periastron. This does not seem physically 
plausible and this evidence, coupled with the X-ray stud- 
ies, causes us to rule out this physical mechanism as an 
explanation for the timing residuals in PSR B1259-63. 
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Figure 2. Pre-fit residuals for the Manchester et al. (1995) tim- 
ing solution. The vertical dashed line indicates the end of the 
Manchester et al. (1995) data set. 
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Figure 3. Post-fit residuals for Model 1, i.e. fitting for pulsar 
position, period, period derivative, dispersion measure, Keplerian 
orbital parameters, and P, P and P to model the timing noise. 
The rms residual is 350 fis. 



2.2 Timing noise 

PSR B1259— 63 is a comparatively young pulsar with a pe- 
riod derivative, P, of 2.3 x 10~ 15 . Such a young pulsar should 
suffer timing noise which is usually dominate d by a cubi c 
term and may contain higher-order derivatives (Lyne 1996). 
For a pulsar with such a period first derivative the expected 
value of P is of order ±2 x 10~ 26 s _1 . (The intrinsic f> aris- 
ing from magnetic dipole spin-down, i.e. assuming a braking 
index of 3, is expected to be only of order 10 -28 s _1 .) Fit- 
ting for P, P and P terms gives the residuals illustrated 
in Fig. 3 and an rms residual of only 350 /is. Parameters 
for this fit are given in Table 1 as Model 1. The value for 
P is within the expected range. Fitting for a P does not 
improve the fit and does not give significant values for this 
parameter. 

The parameters in this fit are uncertain in the sense 
that they are relatively insensitive to adding or subtracting 
an integral number of phase turns at the two periastrons, 
provided that the number of turns added or subtracted was 
the same at each periastron. This ambiguity arises because 
of the large (~ 100 days) gap in the data resulting from the 
pulsar eclipse at each periastron. Quoted parameter-error 
estimates take this uncertainty into account by reflecting 
the parameter range for added or subtracted turns giving 
satisfactory solutions, that is, increasing the rms residual by 
less than 20 per cent. Up to 3 phase turns could be added 
or subtracted at each periastron within this criterion. 

It is difficult to derive an 'activity parameter' (Cordes & 
Helfand 1980) for binary pulsars, especially for long-period 
binaries where some of the timing noise may be taken up 
in the orbital parameters. Having removed the lower order 
terms in the fitting process above, the remaining residuals 
in Fig. 3 are probably the higher order terms of the noise 
process. In particular, the ti ming residuals are very similar 
to those of PSR B1951+32 ( |Foster et al. 1994 , a solitary 
pulsar of similar period and youth to PSR B1259— 63. 



2.3 Secular changes in the orbit (6j, x) and 
classical spin-orbit coupling 

For a pulsar in orbit with a Be star, one may expect changes 
in the longit ude of the peria stron, u, and the inclination of 
the orbit, i, (Lai et al. 1995). The latter manifests itself as 
a change of the projected semi-major axis x = a p sin i. The 
physical cause of these changes is the 'classical spin-orbit 
coupling', i.e. the fact that the spin-induced quadrupole of 
the fast-rotating companion leads to a 1/r 3 term in the grav- 
itational potential which in turn leads to apsidal motion and 
precession of the binary orbit (Kopal 1978, Smarr & Bland- 
ford 1976, Lai et al. 1995, Wex 1998). 

Two unweighted fits for the pulsar-spin parameters, in- 
cluding a P to model the long-term behaviour of the timing 
noise, the Keplerian parameters for orbital motion, and the 
two post-Keplerian parameters, uj and x, are given in Ta- 
ble 1 (Model 2 A and Model 2B). The corresponding post-fit 
residuals are shown in Fig. 4. Again, realistic errors for the 
parameters in the table were obtained by adding (or sub- 
tracting) an integer number of phase turns at each perias- 
tron until the fit deviated from the best fit by about 20 per 
cent. In this case, however, two satisfactory independent fits 
could be obtained by having, in the case of Model 2A, an 
equal number of phase turns at the two periastrons and, in 
the case of Model 2B, one more phase turn at the second 
periastron. 

Let us assume that all the timing noise is modelled by 
P, i.e. that the fitted u> and x result only from changes in 
the orientation of the binary orbit due to classical spin-orbit 
coupling, and focus on the question: can the values for lo 
and x, as given in Table 1, be explained by a reasonable 
quadrupole moment and orientation of the Be star? 

In the case of PSR B1259— 63, the spin angular mo- 
mentum of the pulsar is completely negligible compared to 
the spin angular momentum S of the companion and the 
(total) orbital angular momentum L. The angle 9j between 
the conserved total angular momentum J = L + S and L is 
also conserved (see Fig. 5). The intersection of the orbital 
plane (perpendicular to L) and the invariable plane (perpen- 
dicular to J) has the longitude $ and the longitude of the 
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Figure 4. Post-fit residuals for Model 2A (upper) and Model 2B 
(lower), i.e. fitting for pulsar position, period, period derivative, 
dispersion measure, Keplerian orbital parameters, ui, x and P to 
model the timing noise, using the BTH — (- timing model (see Wex 
1998). The two models differ in the number of turns added at 
periaston (see text). Model 2 A has an rms residual of 340 fis and 
Model 2B an rms residual of 390^s. 



periastron with respect to the invariable plane is denoted 
by "J. The angle between L and the direction of sight Ko 
is denoted by i. The angle i is not uniquely determined by 
timing observations and is equal to either 36°or 180°-36°for 
PSR B1259-63. 

The classical spin-orbit coupling causes a (linear-in- 
time) precession of the angles <!> and ^ which implies a 
rather complicated change of the angles i and u>. For the 
case |S| <C |L| [f] the changes of the longitude of periastron, 
u), and the projected semi-major axis, x, are approximated 
by 

3n 7/3 J| 



" Q ~ 2(GA/) 2 / 3 (l - e 2 ) 2 

«(.-! 

xq_ ^ 3n 7/3 Jj 
x ~ 2(GM) 2 / 3 (1 



sin 9 + cot i sin 8 cos 9 cos $ 



! (1 - e 2 ) 2 



cot i sin 6 cos 9 sin $ , 



(1) 



(2) 



t For PSR B1259-63, |S| ~ 0.1|L|. The approximation \S\ < |L| 
is justified for the accuracy required here. 
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Figure 5. Definition of different angles in the binary system. The 
invariable plane (X-Y) is perpendicular to the total angular mo- 
mentum J = L + S, and the line-of-sight unit vector Ko lies in the 
Y-Z plane. 8j is the inclination of the orbital plane with respect 
to the invariable plane, <3> is the longitude of the ascending node 
with respect to the invariable plane and ^ is the longitude of the 
periastron with respect to the invariable plane. 9 is the inclina- 
tion of the companion spin with respect to the orbital angular 
momentum, ij and i are the angles between the line of sight and 
J and L respectively. 



(Smarr & Blandford 1976, Lai et al. 1995, Wex 1998). n = 
2n/Pt is the orbital frequency of the binary system, G is the 
gravitational constant, and M = m p + m* the total mass of 
the system. The quantity J| is a measure of the quadrupole 
moment of the Be star 

Jl = R-R = \ k R 2 M , (3) 
m* 3 

where fl„ = Q,, j '(Gm, / Rl) 1 ^ 2 is equal to unity if the star is 
rotating near break-up. 7| and I* are the moments of inertia 
about the spin axis and an orthogonal axis, respectively. 
R, is the (equatorial) radius of the Be star and k is the 
apsidal motion constant which is of the order of 0.01 for 
a 10 Mq main-sequence star (Schwarzschild 1958; Claret & 
Gimenez 1992) . Assuming a typical radius of R, ~ 6Rq and 
a rotation of f2* ~ 0.7 for the Be-star companion we find 



J 2 * ~ 0.1i?Q = 2.2 x 10 
The ratio of equations 



6 AU 2 . 
Eh and 



leads to 



UJQX 



sin $ — cos "!> 



1 + 3 cos 29 



(4) 



(5) 



xq 2 cot i sin 29 

an equation that does not contain the uncertain quantity 

ft. 

The observed value of u> (see Table 1) is a compound of 
ujq and the general relativistic contribution 

-i 



Cogr ~ +0.00003 deg yr" 
Thus we have for PSR B1259-63 

( Coq = -0.00021(1) degyr" 
x Q = -0.15(3) x lO" 10 



Model 2A : 



Model 2B : 



6jq = -0.00043(1) degyr" 1 
xq = -2.40(4) x 10" 10 



(6) 



(7) 



(8) 
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Figure 6. 9-& parameter space for Model 2A (upper) and Model 
2B (lower) for the case i = 36° ± 5°. (for i = 144° ± 5° one has 
to replace <I> by <J> — 180°). The grey areas are excluded because 
x < 0. The regions to the left of the left dashed curve and to the 
right of the right dashed curve are excluded because u>o < 0. The 
solid curves border the narrow region of possible values of 9 and 
<&. For the regions to the left of the left dotted curve and to the 
right of the right dotted curve the Be star's rotation exceeds 70% 
of the break-up velocity. 

The restrictions implied by equation (jHJ) are plotted in Fig. 6 
which, in addition, shows restrictions implied by optical ob- 
servations of the projected stellar rotation velocity, v h s : 



± \ 1 — I - — — — ] = cos i cos 8 + sin i sin 8 cos $, (9) 

where v mMC is the break-up velocity. For SS 2883 one finds 
180 km s" 1 for fobs (Johnston et al. 1994). As mentioned in 
Section 1, the break-up velocity for Be stars is not very well 
known; we use u max = 400 km/s as the break-up velocity of 
SS 2883. If we assume that SS 2883 rotates at 70% of its 
break-up velocity (S7* = 0.7) we get the dotted line in Fig. 
6 and allowed values of (8, <E>) are therefore close to (75°, 
355°) and (105°, 175°). 

The large value of 8 implies a significant birth kick to 
the neutron star (cf. Kaspi et al. 1996). Since recent mea- 



surements of pulsar velocities suggest kicks of 200 - 400 km 
s _1 are commonplace (Lyne & Lorimer 1994), there is a high 
probability that PSR B 1259-63 was kicked into a highly in- 
clined orbit; indeed, there is a high probability that an ec- 
centric orbit resulting from a birth kick is retrograde (Hills 
1983). 

Using equations (jjj) and (||) one can calculate J| as a 
function of $, and so using equation one obtains fcP* as 
a function of Assuming that R, ss 6Rq for B2e stars, we 
find those values for the apsidal motion constant k which 
are necessary to explain the parameters u> and x for Models 
2A and 2B by classical spin-orbit coupling exclusively. We 
find that Models 2A and 2B need rather large values for 
the apsidal motion constant k. While the minimum value 
for Model 2A, k = 0.02, is still in the range of theoretical 
models (Claret & Gimenez 1991), the value of k = 0.05, 
necessary for Model 2B, seems to be too large. 

As in the previous section, this model provides a good 
explanation of the data, i.e. this solution maintaines phase 
to better than 3% of the rotational period. The remaining 
residuals are fully consistent with higher order terms of the 
timing noise (c.f. end of Section 2.2). However, it is likely 
in this case that such timing noise cannot fully be removed 
with a P term and that this may affect the measured values 
of ill and x and hence the derived values for 8 and This in- 
troduces an additional uncertainty which is not represented 
by the quoted error limits. 

Simulations suggest that observations around and af- 
ter the next periastron (1997 May 29) should distinguish 
between the various models. 



3 FURTHER EFFECTS ON THE ORBITAL 
MOTION 

Besides quadrupolar effects resulting from the oblateness of 
the Be star, orbital perturbations could arise from tidal ef- 
fects on the Be star, frictional drag as the pulsar passes 
through the circumstellar disk, and mass loss from the Be 
star. We will show that the major influence of these effects 
is to change the orbital period, p,, by at most 3 seconds per 
orbital period (P& $ 3 x 10~ 8 ). However if we additionally 
fit for p, in the models described above, the rms does not 
change (i.e. the value of p, is not significant) and the la er- 
rors on Pb are of order 10 . Hence none of these effects have 
any impact on the current timing solution at a measureable 
level. 

3.1 Tidal effects 

In the case of PSR B1259— 63, perturbations due to tidal 
effects were investigated by several authors using stan- 
dard theory of tidal dissipation and found to be negligible 
(Kochanek 1993, Manchester et al. 1995, Lai et al. 1995). 

Resent work by Lai (1996,1997) and Kumar & Quataert 
(1997) shows that (differential) stellar rotation can change 
the strength of the dynamical tide significantly. In particu- 
lar, retrograde rotation (with respect to the orbital motion 
of the pulsar) increases the energy and angular-momentum 
transfer by more than two orders of magnitude. This effect is 
most likely the explanation for the observed p, of 3 x 10~ 7 
in the main-sequence star binary pulsar PSR J0045-7319 
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(Kaspi et al. 1996). However the periastron distance of PSR 
J0045-7319 from its companion is only 4 stellar radii, while 
PSR B1259-63 passes periastron at a distance of 24 stellar 
radii and thus the strength of the dynamical tides is several 
orders of magnitude smaller. Following calculations in Lai 
(1996) we find for PSR B 1259-63 



-10" 



/50yr\ f T 2 (r,)\ 
\ r d ) \ 0.001 J 



(10) 



where Td is the dissipation time for the dynamical tide and 
r\ — (m*/M) 1 ^ 2 (r pe ri/i?*) 3 ^ 2 where r per i is the periastron 
distance. The function T-2,(j]) depends strongly on the rota- 
tion rate of the companion star and increases with decreas- 
ing values for £l* . Assuming nearly maximum retrograde ro- 
tation of the companion star, Lai (1997) found T 2 = 0.01 
for PSR J0045-7319 (77 = 7). For PSR B1259-63, T 2 should 
be clearly smaller, since r\ — 110 for this binary system; 
one expects T 2 & 0.001. Therefore, for realistic values of 
the dissipation time, Td, the dynamical tides are absolutely 
negligible in the PSR B 1259-63 system. 

3.2 Fractional drag 

A fit to dispersi on-measure changes around the 1994 Jan- 
uary periastron (Johnston et al. 1996) indicates that the ra- 



dial dependence of the electron density in the circumstellar 
disk of SS 2883 has the form 



n e (r) ~ 4.5 x 10 12 (r/7?„)" 



(11) 



Assuming a hydrogen plasma, this corresponds to a mass 
density 



p H {r) ~ 7.5 x 10 _1/ (r/#*) g cm 



(12) 



The radius within w hich m atter is captured by the pulsar is 
given by the Bondi (1952) relation 



2Gm„ 



(wrei /100 km s 



-^2 



(13) 



where v IC \ is the velocity relative to the surrounding medium. 
The drag acceleration induced by the accretion of matter 
onto the neutron star is then given by 



A = - 



(14) 



To estimate if friction has an important effect on the orbital 
motion of PSR B1259— 63, we assume that the circumstellar 
disk is not inclined with respect to the orbital plane, im- 
plying a maximum effect. We further assume (to simplify 
the calculations) that the velocity of the pulsar around pe- 
riastron (~ 150 km s _1 ) is much larger than the velocity of 
the circumstellar material^]. In this case the direction of the 
friction is roughly anti-parallel to the velocity of the pulsar, 
Vrd = v p . Friction will not affect the plane of the orbit, but 
it will affect the Keplerian parameters. For the change of the 
orbital period Pt, eccentricity e, and longitude of periastron 
u one finds (see e.g. Danby 1962) 

t The velocity of the circumstellar matter is expected to be of the 
order of 150 to 300 km s — 1 (Waters 1986; Bjorkman & Cassinclli 
1993). Although we neglect this fact, we will still obtain the cor- 
rect order of magnitude. 



" (2ttGA/) 2 / 3 : 

2 sin tp , 
to = — A , 



. =2 cos^+e 



(15) 



(16) 



(17) 



where tp is the true anomaly of the pulsar and v the abso- 
lute value of the relative velocity between the pulsar and its 
companion. Integration of these equations leads to the fol- 
lowing change of the orbital parameters during one orbital 
revolution: 



APi, « -0.002 s , 

Alo = 0, 

Ae w -3 x 10" 1: 



(18) 
(19) 
(20) 



These values are well below the present measurement pre- 
cision. Thus we conclude, in agreement with Manchester et 
al. (1995), that frictional drag has a negligible effect on the 
TOAs of PSR B1259-63. 

^From equation (p(i[), we can calculate the typical time 
scale for circularisation: 



~ ire ^ -r—Pb ~ 10 13 yr . 



Ae 



(21) 



There are indications that the circumstellar disk is tilted 



with respect to the orbital plane (Johnston et al. 1996). 
Therefore the pulsar is within the circumstellar disk twice 
per orbit, but only for a comparatively short time, which re- 
duces the effect of friction by at least an order of magnitude. 
Furthermore there is no evidence that any accretion of ma- 
terial onto the pulsar took place during the 1994 periastron 
passage (Tavani & Arons 1997). 



3.3 Mass loss of the companion 

If the companion star loses mass at a rate m» so that there 
is no linear momentum loss in the instantaneous rest frame 
of the star, the orbital period will change by a rate given by 
(Jeans 1924, 1925) 



(22) 



P b _ 2m, 
Pb m, + m v 

The mass loss of SS 2883 is of order -5 x 1O~ 8 M yr" 



( Johnston et al. 1996) and so 



■ 3 x 10" 



(23) 



which gives an orbital period change of ~ 3 seconds per 
orbital revolution, a value which is well below measurement 
precision. 



4 RELATIVISTIC EFFECTS 

Because of the high eccentricity and the large mass of the 
companion there are two relativistic effects which should 
be considered: the Eins tein delay described by the post- 
Keplerian parameter 7 (Blandford fc Teukolsky 1976), and 



the Shapiro delay, which depends on the orbital inclination 
i of the bina ry system and is proportio nal to the mass of the 
companion ( Damour fc Deruelle 1986 ). 
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For PSR B1259— 63, the parameter 7 is quite large. We 



find 



7 = 



c 2 \2tvJ 



1/3 G 2/3 m,K + 2m,) 



0.54 s. 



(24) 



(m p + m*) 4 /3 

As pointed out by Blandford & Teukolsky (1975) and Brum- 
berg et al. (1975), 7 can be isolated only in the presence of 
apsidal motion and requires observations over a time inter- 
val in which co changes by a significant amount. Since, for 
PSR B1259— 63, lo changes by a few arcseconds each orbital 
revolution and less than two orbital revolutions have been 
observed, we cannot extract 7 from the present data; it is 
absorbed into the Keplerian parameters x and to. To a first 
approximation we find 



x + 



7 cos LO 
(l-e 2 ) 1 ^ " 

7sinoj 

z(l- e 2)l/2 



: x — 0.84 sec , 



0?031 . 



(25) 
(26) 



Thus the observed values are offset from the true values by 
a few hundred times the errors of the measurements. 

A second potentially important relativistic effect is the 
Shapiro delay, the propagation delay caused by the gravita- 
tion al potential of the compa nion star. The delay is given 
by ( |Damour fc Deruelle 1986] ) 



Ac = 



2Gm* 



In {1 — e cos U — [sin o;(cos U — e) 



+(1-0 



2x1/2 



cos to sm 



U] sin i } , 



(27) 



where U is the eccentric anomaly of the binary orbit. For 
PSR B1259-63 {i ft! 35°) A s has a range of 400 ^s and 
thus is much greater than the error in most of the TOAs. 
The fact that there are no observations around periastron 
where there is a sharp peak in As reduces the actual span 
of As to 150 fis. This remaining effect can be absorbed in 
the other parameters without changing them by more than 
the typical given error (see Table 1). 



Tidal dissipation and frictional drag in the circumstellar 
matter is shown to be negligible. The influence of the mass 
loss of the companion is too small to be detectable unless 
the mass loss is <; lO _5 M0/yr. At the same time we can 
exclude a significant orbital period change in the TOAs of 
PSR B1259-63. 

PSR B1259— 63 should show the largest Einstein delay 
and largest Shapiro delay of all known binary pulsars. The 
small change in the longitude of periastron makes it impos- 
sible to isolate the Einstein delay. The Shapiro delay peaks 
sharply around periastron and is so far unobservable. 

Again, we stress that the physical parameters given here 
for the companion star, 6, $ and k, should be understood 
as one possible explanation for the significant values of d) 
and x in Models 2 A and 2B. If the long-term behaviour of 
the timing noise of PSR B1259— 63 is not fully modelled by 
a cubic term (P), it is possible that rather large fractions 
of these parameters are not explained by a precession of 
the orbital plane but have their origin in unmodelled timing 
noise. The parameters here show that, in principle, all of 
the Co and x can arise from classical spin-orbit coupling, for 
Model 2 A in particular. 

Although Model 1 gives a good fit to the TOAs with- 
out making use of the classical spin-orbit coupling, it seems 
unlikely that the classical spin-orbit coupling is of no impor- 
tance for this system. 
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5 CONCLUSIONS 

At present, the timing observations for the binary pulsar 
PSR B1259— 63 span seven years. Because of the gaps in 
timing observations around the two periastrons and the large 
timing noise present in this young pulsar, we still are not 
able to derive a unique timing model to explain the TOAs. 

Model 1 is a timing solution including a non-precessing 
Keplerian orbit and timing noise represented as a polyno- 
mial of fifth order in time. This model provides a satisfactory 
fit to the data. The remaining timing residuals are under- 
stood as short-term timing noise similiar to that seen in 
observations of other young pulsars (cf. Foster et al. 1994). 

Equally good results were obtained by Model 2 and 3. 
Both timing models contain just a P term to account for the 
long-term behaviour of the timing noise, and Co and x, which 
both are understood to result from a precession of the or- 
bit. This orbital precession can be explained by the classical 
spin-orbit coupling caused by the quadrupolar nature of the 
main-sequence star companion. The corresponding advance 
of periastron is negative and thus the companion should be 
tilted by more than 30° with respect to the orbital plane 
(See Fig. 6). This can be explained by a birth kick for the 
pulsar (cf. Kaspi et al. 1996). 
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Timing models for the long- orbital period binary pulsar PSR B1259- 



Table 1. Parameters for the three timing solutions for PSR B1259— 63, Model 1 (p, p), 
Model 2A (u>, x), and Model 2B (u, x). The data span is MJD 47909 - 50448. To estimate 
the parameter errors we compared all solutions for various combinations of integral numbers 
of phase turns at the two pcriastrons which give an rms residual not worse than 1.2 times the 
quoted rms residual. The error quoted in the table is the maximum difference between the 
values given in the table and the values obtained while changing the phase turns. Compared 
to the errors given by TEMPO and the errors obtained by the bootstrap method, the errors 
given here, in most cases, are more conservative. 



Model 1 



Model 2A 



Model 2B 



Rms residual (/is) 

a (J2000) 
<5 (J2000) 
P (ms) 

p (io- 15 ) 

Epoch (MJD) 
DM (cm~ 3 pc) 

Pb (d) 
x (s) 
uj (deg) 
e 

T (MJD) 



- 1 ) 



P (io- 26 

P (10~ 34 s~ 
"p (IO" 42 s- 3 ) 
ui (deg yr" 1 ) 
x (10~ 10 ) 

Epoch of Third periastron (UT) 



350 

13 h 02 m 47!65(l) 

-63°50'08'.'7(1) 

47.7620537(1) 

2.2785(2) 

48053.440 

146.78(5) 

1236.7238(1) 

1296.4(3) 

138.668(7) 

0.86990(4) 

48124.348(4) 

-5.2(8) 

+7(2) 

-5(2) 

1997/5/29 19 h 



340 

13 h 02 m 47!66(l) 

-63°50'08'/7(l) 

47.7620537(1) 

2.27709(1) 

48053.440 

146.81(5) 

1236.7231(1) 

1296.3(2) 

138.667(6) 

0.86989(3) 

48124.350(4) 

-0.884(9) 



-0.000184(8) 
-0.15(3) 
1997/5/29 19 h 



390 

13 h 02 m 47!65(l) 

-63°50'08'.'7(1) 

47.7620537(1) 

2.27720(1) 

48053.440 

146.85(5) 

1236.7244(1) 

1296.4(3) 

138.670(7) 

0.86991(4) 

48124.348(4) 

-0.788(1) 



-0.000396(9) 
-2.40(4) 
1997/5/29 19 h 



